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ABSTRACT 

A satellite whose orbital rate is commensurate (or rational) with the 
earth’s rotation rate will describe a constant ground track over a whole 
number of sidereal days. On such a satellite, the small orbit averaged 
along track force from certain longitude harmonics of gravity can build up 
over many commensurate ( synodic) periods to produce appreciable changes in 
the semimajor axis of the orbit and longitude placemeht of the ground track. 
The libratory nat-ure of these changes and the specific harmonics which can 
cause them for any resonance case are discussed with reference to a circu- 
lar orbit. General formulas are presented for these librations, with typi- 
cal errors determined by numerical integration. Tlie specific theory of the 
24 hour, nearly circular orbit, satellite is presented in detail with actual 
orbit data from three such satellites used to determine the resonant earth 
gravity sectorial harmonics through third order. These are found to be: 

It is shown that a number of subsynchronous circular resonant orbits 
offer similar promise in discriminating easily and uniquely other low and 
high order longitude harmonics; among them H54, 
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INTRODUCTION 


It has been known for some time that satellites whose periods are 
commensurate with the earth’s rotation period can suffer long term per- 
turbations in their elements due to longitude variations in the earth’s 
gravity field [l]. The reason this is so is that viewed in a system of 
coordinates rotating with the earthy the satellite with commensurate 
period maintains a fixed relationship to the underlying longitude gravity 
field (Figure l) , Assume there is net change in an element over the 
period of commensurability, due to a gravity perturbation* Then^ no 
matter how small the periodic perturbation is^ over many commensurate or 
synodic periods a large change in that element can occur as long as the 
geographic configuration of the orbit is relatively unchanged. The 
large^ or secular^ change is merely the sum of the many small net 
periodic changes. In this paper I will be concerned only with developing 
the secular effects of longitude gravity harmonics on the semimajor axis 
of circular orbits with constant or nearly constant ground track. The 
orbit period of such a satellite is said to be commensurate or resonant 
with the earth’s rotation period. In the circular orbit case^ without 
inclination restriction^ it is possible to derive the secular change of 
the semimajor axis from a very primitive viewpoint by a straightforward 
orbit averaging of the perturbation forces. A simple investigation of the 
symmetries of the harmonics with respect to the constant ground track will 
be shown to be almost .sufficient to discriminate those. earth perturbations which 
can cause secular changes of the semimajor axis and those which cannot as 
long as the orbit remains circular. As might be inferred^ the long term change 
of the semimajor axis with respect to the longitude harmonics of gravity 
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is strongly dependent on the longitude orientation of the grovind track. 
This may be specified, for example, by one of the equator crossings of the 
satellite. Since a change in semi -major axis is accompanied by a change 
in orbit period, the longitude orientation of the ground track will also 
change in response to the long term perturbation. The residual pertur- 
bations themselves are strongly dependent on the orientation of the 
ground track with respect to each longitude hamonic. Thus, we may ex- 
pect the evolution of the orbit under resonant earth gravity to be 
essentially described by a system of two coupled equations in the longi- 
tude orientation and semi -major axis of the satellite's orbit. We will 
easily derive close approximations to these equations for any nearly 
resonant orbit. One of the critical questions in the use of these equa- 
tions is how many relevant earth harmonics must be Included in any appli- 
cation to get a good representation of the long term evolution. Where 
no knowledge of the longitude harmonics exists, an infinite number of 
them must be carried. But in the case of the earth, it is shown that for 
the high altitude 24 hour satellite (with period equal to the earth's 
rotation period) it is highly probable that only resonant harmonics of 
second and third order need be retained to maintain acceptable accirracy 
in the equations over long periods of time. 

As a reverse application of the resonant orbit evolution equations, 
we solve for the underlying resonant longitude gravity field of the earth 
from actual position data on three synchronous satellites over a three 
year period. Further investigation of these equations using an earth 
model based on recent data reveals that a number of other resonant orbits 
of less than 24 hours may show sufficiently strong long term effects to be 
useful in discriminating precisely other longitude gravity haimionics. 
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THE GROUND TRACK OF RESONANT ORBITS 




Figure 1 shows the groxmd track for a general resonant orbit of 
moderate inclination. 


GltEE/OWiCH 





Figure 1. Gromd track of a Resonant Circular Orbit 

The resonant satellite makes f orbits ( 0,^) over radians of 

geographic longitude (with r and ^ whole numbers) before the gro\md track 
begins to repeat at, say, an ascending equator crossing of from Greenwich. Now 
consider the geometry of the orbit in inertial space (Figure 2) , At time 
t, the satellite is at S , from its ascending node. The inertial longi- 
tude from its node, AL , is given by 
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Figure 2. Orbital Geometry of Resonant Satellite in Inertial Space 


Let us assume a prograde orbit, with the earth's orbital rate. 

Thus, the geographic longitude at S is: 


X= Ao -HAL 


( 1 ) 


In the resonant or constant gromd track orbitjwhen X»A© — AX s= 

. Let be the satellite's orbital rate, and 

y where and are arbitrary real n-umbers 

at this point in the analysis. Also, let zl ^ the longitude 

turned by the earth since time zero. Between ground track repetions, 

I I 

with M in days. Then^ since: 


(la) 


(l) becomes; 


A A = ^ -i) 


( 2 ) 


E^>] 


in (1") is zero and periodic every zr;Lir'r- *^4 


Therefore, the ground track repetion condition: 


A\- j 


, implies 
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From ( 2) : 


(at repetition of ground track) 

= = A7Tn’(^ -1^ 




or 


( 3 ) 


Also since _ ;2.irr at ground track repetition. 


i = <*4 

«• r “e ^ 

which means that and must be whole numbers for resonant 

orbits. With a given ratio, n* is the resonant, repetition, 

or synodic period in sidereal days and I" is the number of orbits in the 
resonant period. Eq,uation (4) in ( 5 ) gives 

» 


t 


ss r-H 


ih) 


( 5 ) 


Thus, from ( 5 ) and (4), a specification of and giving the 

commensurate orbital frequency^ determines T from n’ as a least 
common denominator of (4), and then ^ , the global circuits in the 
resonant period. If ir>H* , fx^ft satellite period 

is less then 24 hours. If 
satellite period is greater than 24 hours. 
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RESONAIW ORBIT EVOLUTION 


Consider the perturbation forces on the resonant circular orbit. 

The only component of these forces which can effect a change in the semi- 
major axis of this orbit is the along track component Fy (Figure 3) . 

This is so because the instantaneous semimajor axis of an orbit is only 
a function of its total energy^ kinetic and two-body potential. The only 
force which does work on and thus effects the total energy of the circular 
orbit satellite^ is the along track force. 


So 


Figure 3. Along track perturbation force on a circular resonant orbit 

In an elliptic 
satellite is: 


where CL is the semimajor axis. To first order^ then^ the change in the 
semimajor axis due to a small change of energy AH is^ 


orbit of the earth the total two body energy of a 




( 6 ) 
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From (6): 


A<J- = . 


( 7 ) 


But for a circular orbit satellite, the work done by a small^ orbit, 
varying, along track force adds energy the satellite 

each r orbits, amounting to: 


(p er r orbits) S iFy j 


( 8 ) 


where fv is the orbit averaged force defined by: 

•r u t 

' T - 2 


( 9 ) 


Combining Equations (8) and (7)^ the evolution of the semimajor axis of 
a nearly circular orbit under a small perturbing force, must be governed 
by the difference equation: 

t F 

-A A. (per r orbits) = 'T* 

In the resonant case, t~ orbits are made in n‘ sidereal days so that in 
time units AT of n' sidereal days (the synodic period): 




(lo) 
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Considering 4 A- small in terms of the long term orbit 

evolution^ (lO) can be vritten as a differential equation with respect 
to time in sidereal days ( -flX — 


i = 47Tafr ^ 


length mlts 
sidereal days 


(11) 


It is to be understood, of course^ that equation (ll) is an orbit 

averaged differential equation. To derive the orbit averaged longitude 
motion of the initially constant ground track under a sustaining orbit 
expanding force , we resort to Kepler’s period law: 




■^k 


(lla) 


Thus^ to first order^ if the semimajor axis changes by the period 

changes by 

But (lO) in (12) gives the period change due to fy after one synodic 
day as : 


AT (in one synodic day) 



If, in Y~ orbits, the satellite’s period changes by ^ ^ where 


(13) 
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is in units of sidereal days, then the satellite retums to its ascending 

equator crossing in days. The earth meanwhile has 

turned 2 .-TT( rad/day ) [lri%-rAT 2 ] (days) radians. 

But since \r }is K1 , a whole number, the equator crossing is 
zrrrAij, west of where it was at the beginning of the resonant 

period. Thus the drift rate of a resonant orbit after 1 synodic period is 


AX 


- 2.TT r A'Q 


But • 


= 


AX ( any units ) x Tj — AX iL 

T T r 


Thus 

j 


^ ^ (at one synodic day from resonance) — 

AX 


^z.irv\ AX 

r 


Substituting (lla) and (13) in (ih) gives the drift rate with respect to 
the synodic period at the first synodic day as s 


A _ _ iz.TrVn‘/v 


rad/synodic day 


(ik) 


(15) 


In the second and succeeding synodic days^ if the drift rate and pertur- 
bations are small^ the continually acting residual force will build up 
the drift rate linearly with respect to the synodic period^ according to 
(15) whose elements may be treated as constants. Thus (I5) also repre- 
sents^ approximately^ the resonant^ long terrp acceleration of the geogra- 
phic longitude of the initially constant or nearly constant ground track of 
the satellite. Or^ treated as a differential equation with respect to the 
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elements at successive synodic periods 


** 7 * P" 

X 3: ^ IZTT m Ft 

But n* sidereal days = 1 synodic day. 
Thus (16) can be rewritten: 


radians 
synodic day^ 


(16) 


X = Ft 

The relevant (longitude and latitude) earth perturbation forces (per 
unit mass) can be conventionally written as F = F 

where Ck. is the initial or instantaneous radius of the circular 

• • 

resonant orbit. Writing simplify (ll) and (17) 

to: 


radians 

2 

sidereal day 


( 17 ) 


= 4TT ( ►-/«■) 


l/sid day 


( 18 ) 


X 


-l2.TT^(»7n') Ft 


rad. / Sid. day^ 


(19) 
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EVALUATION OF THE RiSONANT PERTURBATION FORCE IN TEM^ OF GROUND TRACK 
COORDINATES 

We can now proceed to evaluate the orbit averaged perturbation 
forces. The simplest way of doing this appears to be in terms of the 
ground track coordinates latitude and longitude which are the coordinates 
of the conventional spherical harmonic expansion of gravity given below. 

We take the gravity potential on the circular orbit as the conven- 
tional associated legendre expansion: 

Ve, = / I - Z. II ( Pn f ^ ) 3;„Cosi*i(X-X„ 

n=l W-6 J f (20) 


where Rq is the mean equatorial radius of the earth, ^ is the latitude 
of the satellite and A is the geographic longitude. By convention; 


. The associated legendre polynomials 


are defined as follows [2.]: 


p% = 0-x^) J- 

» M — 

JLtT 


J 


where 


(20a) 












(20b) 


A. Zonal Gravity Forces 

The zonal gravity potential ( ) is not dependent on longitude. 

The latitude dependence of this part of the potential can be expressed as; 
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( 21 ) 


With the series terminating at the term in if n is odd, and in 

a constant term if n is even# Consider a typical single orbit 

ground trace of the resonance case in Figure 1. 
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Figure h. Orbital Ground Track vith Zonal Perturbation Forces 


The latitude force arising from (2l) is 

a 

tude dependence is: 


Its lati- 


‘ ^ ^ f 


t C 9 ^ 

P =z co^S )a. -f- a. 45/A»<f) 

'n,oflo ^ L ' 

r / M-/ , W-3T ^ 


(22a) 


(22h) 
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9 

‘f >1 


is an odd function and 


,000 is an even func- 


tion with respect to the equator. With respect to the longitude of the 

descending equator crossing of the track, the azimuth t>C is an even 

/ 

function, but F* - is "odd" (See fig. 4). Thus the in- 


cremental work done 




lon the satellite 




) is an odd func- 


tion with respect to the longitude of the descending equator crossing. 


The integral of AW ^ or the total work done by 




over an 


orbit must therefore be zero^ as it must be over orbits^ or the 
"resonance" or synodic period . As long as the orbit remains circular 
the even zonals can have no long term influence on the longitude accelera- 
tion of the constant or nearly constant ground track. On the other hand, 
for the odd zonals, we can use the property of the ground track that 
cos (azimuth) is antiperiodic over , or half the orbital track. 

Since Fh^ood is periodic over 5 ^ 

— COO COS (azimuth)^ is antiperiodic over the 

incremental work due to the odd zonals is also antiperiodic over 
or the total work due to them is zero over any number of single orbits. 

Concluding^ as long as the orbit remains circular, the perturbing 
zonal harmonics of gravity can have no long term influence on the longi- 
tude acceleration of the resonant orbit. 
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B. Longitude Gravity Forces 

With respect to the periodic longitude gravity forces, since cosfl( 
is antiperiodic and sin«( is periodic over a half orbital track, we voiild 
like to find those longitude harmonics whose latitude force components 
are antiperiodic and whose longitude force components are periodic over 
, With such forces^ the long track component is period every 
half orbit^ orbit, and so forth to r orbits or the synodic period 
(See Figure 5). There does not seem to be a possibility of the along 
track force having periodicities or antiperiodicities over spans less 
than Xy^ throughout the resonant ground track (but see Appendix A), 
Therefore, in the general case for an arbitrarily inclined orbit and 
longitude placement it appears that the mean value of these periodic 

fluctuations (when found) will not be zero. In fact every detailed cal- 
culation of the general case of a specific longitude harmonic for a 
resonant orbit has shown so far that only when there is periodicity in 

over , can the orbit average of tie different from 

zero. An example of this calculation, for the orbit averaged effect 
of the H32 harmonic on the 12 hour resonant orbit ( = M sf ), 

is found in Appendix C. These detailed calculations suggest that 
can always be decomposed into the product of orthogonal functions in the 
intervals or ^Xyj^(a full orbit). Such a decomposition would 

allow us to make firm, instead of merely suggestive conclusions as to 
the relevancy of any particular harmonic on a particular resonant orbit. 

r 

However, in those cases where we can show antip eriodicity of \ y^ 

or ^ it is evident that those particular harmonics are not 

relevant to the orbit evolution we are interested in. 
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Figure 5 illustrates the method we will use to show periodicity or 
antiperiodicity of over \,h. In achieving periodicity of 

i ^ 

over (from to Q) in Figure 5 ), we can 

follow the changes of the longitude harmonic forces first from ©tod), 
a constant latitude path, then from© to 0, a constant longitude path. 
This scheme is suggested by the form of the longitude and latitude force 
components of the longitude harmonic . It may be readily 

shown (following the development in Ref. [i]) that the longitude and 
latitude dependencies of these force components (derived from the 
gradient of Eq. (20)) are: 


F' 


+ £05“^' Sl/Jf^b’si/J f f t-bi j COSix(X-AnM)^ 

( 23 a) 


and 


KC4>A) = 

L. ^ c J 

Poh y^-M AViei J 




( 23 b) 


F CifA) 


^ -f <f> -t- • • • CoS ^ gftJ ^ ) 


(2l<.a) 


AtJO 


- 16 - 






Co5***<f <* ^ ■** » • • t 

(24b) 


For n-m odd. 

^lovy, CoSw\CX-XnM^ a.nX StiJt*Lh.^X<^wi) A,re penoilc oVer X*/a, C *•'*'/ 
sCosM(X-A^.-,\i only i a«i 5•^» =o . -rUese 

(for m o.rt 5*.^» 5 f»c.«A ovil^ ^or \t/^^o (a,ll peii-tiv^ or X*/j, 
a.V“bi'fc^'ft.ry A.nAr 


vr» i jni. f tx,±^^ i4 . . n , 

X./a."* 


(25a) 


5ii«^i )tt.rl^^ Cos i^ o.n<i (,X-Ai»iimv^ tftre Awti par’»o<l'»e. Co5i*»(XtX»^jj”Xi*H^ 

^ -cos^^^CX-X^^M^ cos iwX*/ji,= -/ ft-ndl SWMX.yf^ss <?, "rUese. 

coM<I.Tii&vis o-re 5 <ci’*s^re«^- X ‘/g, • 


Wv ■=: »' 


X»/( 


r*!- . 


a. 


Further, examination of (25) shows that Z4 w-wt Is cVen : 


(25b) 






and 


'Aywm A,ni** 

Similarly, (24) shows that if n-m is odd 

PiCt) = Fj, . 


(26a) 

(26b) 


(27a) 
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and 


(27b) 


< 

Thus^ if n-m is even, and we want periodicity of ^ mm over 

(26) shows we can only achieve state(^froin(^by way of the a) and c) forces. 

But since the a^ and c) forces can only be achieved from state if (25a) 

I 

holds, we have the following general result: periodicity of Fy ^^^^over 

a half orbit is achieved for n-m even, if and only if or , 

^'/SL 

Similarly, if n-m is odd, (27) shows we can only achieve the periodic 
state ^ from by way of the b)and deforces, which can be achieved from 
(^only if (25b) holds. Thus we have the alternate result that if n-m is 

I 

odd periodicity of E over Xin is assured if and only if 

. But from( 3 )an' 3 . 0 t]|> bhe half orbit excursion *5; 

HiL ' " ' 


( U -i') = tt-Cj-m/,-). 

Thus the strong presumption is that for a satellite which repeats its 
ground track every r circular orbits over n’ sidereal days, the gravity 
harmonics Hm»v\ which are capable of producing a long term change in 

the period of the satellite are given by: 

V) - M -fw*" _ _ _ 

\^/v ~ ^ ^ ^o\r V1-M &\f&n cKwJi t«.U 

C I — M*/r) Lnifctr*.l PoS'iti^fC >vi 

U'izir) 


(28) 


(29a) 


(29b) 


A 


y\JL 




Ltij fi, - ^ 

Ci-y^*/r) pos'iiwe- KA 

Cvt*=f^ tr ) . 


(29c) 
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Condition (29a) is just that for the 2 k hour satellite (see equation ( 4 )). 
In fact^ for this satellite vith a simple reversal of the periodicity 
arguments above ve can show rigorously that all the remaining longitude 
harmonics, oJ-tL , can have no long term effects on the period 

of this orbit. Similarly, to guarantee antiperiodicity of Fr over 
for any resonant circular orbit, (enabling us to ignore the influence of 
for these purposes) the necessary and sufficient conditions 

are that: 



lie -for 


(50a) 


itX ±4 ) 

!. \n~fA cdtL c,Ael 6>.ll 



C 1 - w'/r) 


(30h) 


C t t iO-*' 

C 1 

^ -for y>-M e,Ve-tn ckaJL ctU ^ 
|*csn 1**^ (.1^ 

Ooc) 


We can go at least one step further in eliminating additional harmonics 
from consideration by utilizing the periodicities of both cosoC and sinoC 
over a full orbit Since a full orbit returns the ground track 

to the same latitude, we only need examine the antiperiodicities of the 
longitude dependent factors of the harmonic forces, cos rv» and 

sin M U-x„ . Rewriting (25b), we see such antiperiodicities are 

possible for 2 . X. ^arbitrary ^ only if: 


since 

s C t M from (28), we have the general result that 
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additional "non resonant" harmonics on a constant groTond track orbit 
(in the sense considered here) are those Hhka for which: 


>v\ 





In passing we can note that ( 3 l)shoiild also be a valid discriminant of 
non resonant hannonics on a constant groimd track elliptical orbit as 
well. This is because (assuming negligible change of all elements in a 
single orbit) the flight path angle is still periodic over 
The tangential component of the radial gravity pert-urbation forces [all 
with longitude dependence cos ] are thus antiperiodic in 

S-Xy for all m given by ( 5I) . 

The criteria (30) and (51) enable us to ignore a large class of 
harmonics on the resonant orbit (specified by with the lowest 

common denominator) . But detailed calculation shows they appear to ex- 
haust the ignorable harmonics for only the 12 and 2 h and 36 hour orbits. 
It can be verified that a single simple formula equivalent to the pre- 
sumably resonant harmonic criteria ( ^ , at least for 

(synchronous and sub synchronous orbits) is: 




= Jl. ( n ^ : p =. o, X. t * • 

n' 


where m and n are positive numbers and p is such that II M 
(see appendix A for a discussion of the supersynchronous ^ ^ 

resonance selector)^ 8ut (32) is Just the resonance condition given by 
Allan [3] for subsynchronous satellites of daily ground track repeti- 
tion ^ r: 1) . Allan *s theory (in terms of Lagrange's planetary 


equations^ orbit averaged over a resonant period of a day) thus appears 
to be readily extendable for a resonant orbit of any 

commensurate period by way of (32). (see Section 5). Previous to Allan’s 
development the author had calculated 'throu^ for 

the synchronous satellite directly from the spherical harmonic potential 

C4-1 

(20) . . The results vere (for the n-m even, relevant harmonics); 

J ^ 


— r r 




SI 




31 




'Ct+coil)lx ^ 


e 


31 




f; 


^ r= — 7 C H-t-osZ 

j 




4 - 


(vhere = semimajor axis in earth radii) . 

Utilizing the same technique of direct trigonometric reduction for 
T-Tw. on the 12 hour satellite , the 

author has shown (see Appendix C) that through 4 th order (n= 4 ), the only 
relevant harmonics are H32 and H44 with: 


c‘ (.12. = ~l£'y^^StlJl(l^'XCoil^3coS^l)coSa.C\-Xiz) 

■T;3a. s — a? 




(53) 


( 54 a) 


22 - 


Similarly, for the 8 hoiir satellite through Jrd order (n= 3 ), the only 
relevant harmonic vas found to he H33 with: 

~ 4^ — 3^ C 1+ 6osZ)5//^3( K-X33^ • (35) 

S 
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EVALUATION OF THE GENERAL RESONANT, ORBIT AVERAGED, ALONG TRACK FORCE 
IN TERiyiS OF KEFLERIAN COORDINATES 


In the previous section, though we developed a niimber of criteria 
hy which we could reject a class of nonresonant harmonics, we could 
not accomplish the two necessary and sufficient tasks for a truly suc- 
cessful resonance theory: 

1. To find the single necessary and sufficient criteria which 
determines to be either resonant or nonresonant on a 

given commensurate orbit, 

2, To determine efficiently by a simply evaluated formula, the 
inclination function or integrals representing the orbit 
averaged behavior of the resonant satellite. 

We will see in this section that through the expression of the pertur- 
bation force in terms of Keplerian elements, instead of ground track 
elements^ these major tasks can be efficiently accomplished. 

We have seen from the previous development (Equation l8) that the 
nondimensional semimajor axis ( ) of the resonant orbit evolves at 

a rate, in units of l/sidereal days, given by: 

d„ = A-lTL >-/„■) Ft ^ 

Where Fy is the K* orbit averaged along track force due to a dis- 
turbing harmonic as a fraction of the principal, radial, gravitational 
attraction. Over one synodic period (n* days), the evolution of the 
semimajor axis is thus given by: 


1 = 4 TT ^ r F-f . 


semimajor axis units 
synodic days 


( 36 ) 


But the instantaneous evolution of the semimajor 

C^XCsJ 

axis is given also from Lagrange’s planetary equations^as: 


- 2k - 


Jia. 

2 t 


ix* 


(57) 


In Eq.uation (57)^ is any disturbing potential function, whose 

\/ ^ 

gradient, by definition_,is the dis.t^lrbing force. ^ is the mean 
anomaly at t=0 (epoch) (or the modified mean anomaly), defined from the 
mean anomaly M by: 

t 

M = + X‘ 

''o 


( 38 ) 


Since the orbit we are considering is circular, we can arbitrarily 
put W (the argiment of perigee) equal to zero. Then M is also the 
true anomaly counted from the ascending node^ or 0 (see Figure 6). 

Thus^ 

t 

^ ( circular resonant orbit) = + x‘ 

o 

We note in passing that from its definition above^ for the resonant 
circular orbit 

^ ( 39 ) 

At any time in the dynamics, all other elements being held fixed. 

Separating ( 57 ) and integrating over I' orbits in n' days, ;mder the 
ass\unption that <V*i a. ) changes negligibly in this tlme_^we have 

Aa.(in n* IamsN « 

^ IX* ^ ^^0) 
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y 


■t^€. $ ti^r'lskil^ po ir&in ^ieul a 

But from Kepler’s period law; 



orV)i. t 


<xvero.^e<i ^erivA.irit^e 


(41) 

(42) 


With the time units of in sidereal days. Equations (4l) and (42) in 
(4o) gives over a resonant time period OuSZ 




semimajor axis units 
synodic days 


(43) 


Equating (43) with ( 36 ), the nondimensional \T orbit averaged along 
track force is given from the averaged potential derivative, by; 


F* = ,.£=.1^ . 


(44) 


Now, in order to find , we have to express the disttirbing 

potential^ in our case^ that arising from a longitude gravity harmonic ^ 
in terms of a set of mutually exclusive orbitial elements including \ 

C33 


Allan has summarized previous developments [ 6 ], [7] of the longi- 
tude harmonics into series in terms of the ordinary Keplerian elements 
OLj 4 .^ W and M or f (the true anomaly). Here ve only 

state the result for zero eccentricity. Each longitude spherical harmonic 
C»^ =/’•<>) in the expansion ( 20 )^ gives rise to a potential func- 
tion^ in terms of the elements of a circular orbit (see Fig. 6 ). 
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rvi (^i/i — **‘^nrvT)^ 

where J — 'J^ , R i signifies the real part °f | j 

and the dynamics (epoch) begins at zero sidereal time, or when the 
Greenwich Meridian passes through the Vernal Equinox. 


(45) 



op Hum 6MV»Vr HAPMcPtC. 


Figure 6. Circular Orbital Elements at time t 


The inclination functions 


F 40 

' MMp 


are given as: 


F = 

* n»\»p 


• S)/J n/a.') 


do 5 itAc). 


( 46 ) 


where 

Wow since 


for the circular orbit. From (45): 


- 27 - 







f= 


=t> 


+* wi ^t/L ^ X . 


^ 47 ) 


We can include the imaginary part of J f~HMp in the exp, term 

of ( 47 ) by defining a phase angle such that: 


» . W-t»^ 


J-J 


— J 


W-M+-I 


— e. 


>AtM 


■“ C.05 ^Mma "t" J 


5iaJ S 


Am 


(48a) 


We see from its definition in (48a) that: 




31 n-M + l^ • 


(48b) 


With the result of (48), the real part of (47) is readily evaluated 
and the potential derivative becomes 


2>M 


— -C^ ) (n-2.p') 

(49) ^ 


VI 


^ X 


where 


F il') — ^VIwp/jA-" 

MiAp ' 


( 49 a) 


Taking the average of over orbits in n* days: 

^x* , 
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assuming constant A. over the resonant period. Let be the 

nodal argument at tsO. Then 6 = , assuming a constant 


orbital rate over the resonant period. In terms of , the 

n+1 integrals of ( 50 ) are: 

'Jit + 

( 51 ) 

Let ^ Vw (t/l - -pShm ^ cKtul (52a) 

( 55 b) 


Under the basic assumption of constant orbital elements over the 
relatively short resonant period n*, and %'xp will be 

constants and the integrals of ( 51 ) can be evaluated as: 



But ^ and are simply related by the constant gro^ind track 

condition^ Equation (4-): 


= ’'A ■ 


(54a) 
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OJs 


with in units of radians /day. 

From (5^a), (54b), and (52b), 

= ^TT - mm'J 

sE joric' 


where 


= rin-a.p) .««' ^ 


o tl, 


depending on the integral values of n, m, p, r and n* for the orbit 
and harmonic in question. 

If n’Skp = ilTr( ±1, ±3^^...^^and ( S,p t i i/J 

In this case then all the n+1 integrals of (50) are zero . However, 


I < ( ^ 

when K - h =0 


and each integral solution from (55) 


is an indeterminate form. To resolve this case we seek formally: 

Wmlt: *777 ^ ^ ^ 

fe-0 ^ ^ ^ ^ I I 

i:.m; t — "3 ia.it ) < , 

By L*Hopital’s rule, this last limit is the limit of the ratio of the 
derivative of the numerator to the derivative of the denominator with 


respect to 
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for can produce residual or long term orbit 

averaged along track forces. 

But this is precisely the resonant harmonic condition presumed in 
(32) and here shown rigorously. We have thus shown that Allan's one 
day resonance theory [j] ^ for circular orbits is readily 

extendable to any number of synodic days. (Tlie extension is made by sub- 
stituting ►"/•I* for p in Allan's results). 

To complete the calculation of the averaged along track force for 
any constant ground track orbit , is resonant 
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on it if a single positive (or zero) p exists which satisfies (58). 

Given n, m, r and n', of course, there can he only one p which does 

satisfy ( 58) . Thus only one of the n+1 integrals of ( 50) can be non- 
zero and with this single resonant n, m, p combination [satisfying (58)] 
in mind, the r* orbit averaged disturbing function due to (resonant) 

is [with the results of (57)^ and (58) in the form n~lf> s j 

(53) and then (50)]: 






Finally we would like to express this disturbing function in terms of 
the longitude of the constant eq.uator crossings in the resonant 

ground track, (see Fig. 7). 



Figure 7, Orbital Parameters at t=0. 


The satellite was at its ascending equator crossing just previous 

_• i / 

to the epoch at S when the Greenwich Meridian was at ^ 
west of Y . We assume, as before^ that the orbital parameters 
(with the modified mean anomaly) suffer negligible change over the 
resonance period. Then since ~ j 

Thus X , the longitude of the ascending equator crossing just 
previous to the epoch is given by: 

+• ' (6o) 

Equation ( in (59) gives the averaged dist-urbing function as: 


>X* 




(61) 


Finally^ it makes no physical difference in our arguments whether we 
reference the orbit to any of the ^ distict^ constant^ ascending or 
descending^ equator crossings ( ir longitudes in total). The K" orbit 
averaged disturbing force on the constant ground track is obviously 
independent of the specific longitude or orbit time reference. 

We can find the representation of the force with respect to the 
other 7}T-k crossings by finding the relationship of these crossings to 
A ^ the ascending equator crossing just prior to the epoch. 
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From (5) and ( 4 ), the constant ground track extends — ^ 

=: 2 .Tr(r-n‘) radians in geographic longitude. Thus^ since there 
are equally spaced equator crossings in this span, they will be 

separated by (, 1^— W } radians. Thus the successive equator 

crossings^ from the original one calculated in ( 6 o)^are give it by 

= X ^ 

Where >C s (odd for descending and even for ascending 

crossings), Solving (62) for X and substituting this in the cos. 
part of (6l)^this part becomes: 

— ^ wittX C ^ ” *^/i" ^ * 

i 

But from the resonance condition (58), s y\~xp , which 

must be a positive Integer. Thus Vv| TT^ ( I (.»*♦' M ^xp), 
which must, for the resonant orbit, be some integral (negative or positive) 
multiple of TT , 

We see from this that for X even, or the ascending equator 
crossings^ ^<>5 <*17C^ the force is expressed 
exactly as in (6l). Furthermore, if { | is even (or zero) the 

force is also given by (6l) with respect to all the descenaing equator 
crossings. Only when ) is odd will the form for the 


% 


force with respect to the descending equator crossings be different 
than (6l). Since ZP is always even^ and ^ this latter 

condition is fulfilled if n-m is odd. We can cover all cases by the 
form (6l) if we write a new phase angle which varies 

not only with C >1/ but with whether the descending or ascending 
node X is being considered. 

Let 




( 65 ) 


Then we can write ( 6 I) in the general form: 






(6k) 


where , I refers to the geographic longitude of any descending equa- 
tor crossing^ and ,X JS 2 - refers to any ascending equator crossing 
of the resonant orbit. Summarizing these results, from (64) and (44), 
the nondimensional 1 “ orbit averaged along track force on the resonant 
orbit tH'r) is given as the sume of all relevant n, m terms discrimi- 
nated by ( 58 ) : 

rele.Vtt.rtt (5^) 
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Where: 




*Sfe).5r'( . )UJ“> •« >' “« • 


ic s H-W , 

MAA 


_ M Crete VA.nt'^ - wt»~e'eVA.r}t^( ►»/►-) 

'' = ^/l^y*** ^ (65c) 


there being only one integral or zero value of ^ [from (65c)] 
for each relevant H»)m ^ resonant circular orbit^ €K.nX 1 


JL^ I 
/si 


for descending equator crossings^ 
for ascending equator crossings^ and 




(65d) 
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ORBIT EVOLUTION EQUATIONS FOR CIRCULAR EARTH RESONANT ORBITS 


Using the results of the last section on the IT orbit averaged 
forces due to the longitude harmonics, the long term semimajor axis 
change of the resonant orbit is given from (65) in (18) as: 

C tt Uiii.tX)/5i^eretcl iaij • 

Siinilarly^ the T orbit averaged geographic longitude C X/) drift 
resulting from the resonant forces is given from (65) in (19) as: 

^ rdeVant ^ 

The parameters and are defined in (65a-d)^ as is the 

inclination factor • 

It should be noted that ( 66 ) and (67) apply strictly^ oniy to an always 
circular orbit satellite. Tliis condition is, of course, immediately 
violated as soon as the perturbations which drive these equations (as 
well as other perturbations) are introduced. In addition, the equations 
were derived with reference specifically to the perturbations from 
the condition of constant ground track. They may be thought of as 
giving the leading term in an expansion for the drift about the reso- 
nance condition. With respect to orbits whose equator crossings have 



an initial drift rate^ these equations, as stated^ are in error even 
at time zero. However, niimerical trajectories have confirmed their 
adequacy in predicting the orbit evolution of many near resonant 
satellites with appreciable drift rate and moderately small eccentricities [8], 
(Allan's paper [5] gives more information bn the magnitude of the orbit 
acceleration arising from nonzero eccentricty, nonzero drift rate, the 
oblateness of the earth, as well as sun and moon accelerations.) The 
most important feature of these equations is that, considered either as 
evolving from a constant orbit or near constant orbit state, only the 
variables time Lofi X* 3 and X change appreciably from any initial 

values. This can only be appreciated by integrating ( 66 ) and ( 67 ) 
from initial conditions under this assumption and checking the change 
in the other variables (principally CL ) of the equations. This has 
been done in the case of the 24 hour satellite under H 22 perturbations 
[9]rTKese tests show ( 67 ) to be essentially a nonlinear equation in X 
with constant coefficients determinable from initial conditions. The 
coupled equation in the nondimensional semimajor axis change ^ 

also may be thought of as having essentially constant coefficients. 

With these remarks, (67) can be reduced to a simple pendulum equation 
[for Xtt)] for each relevant harmonic Depending on 

the initial conditions of the near constant ground track, with respect 
to each relevant harmonic, the possibility of either circulatory or 
libratory drift of X exists, coupled with long period oscillations of 
the semimajor axis described through (66). The characteristics of this 
*long term resonance" or libratory behavior is discussed in Section 9- 
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It has also been discussed in the literature previously in references 
[lO]^ [ill and ri 2 ] for synchronous satellites and Ref, [13] for 
other resonant orbits. Here^ we would just like to estimate the 
maximum long term effects the longitude harmonics may have on near 
circular subsynchronous resonant orbits which may have futiure geo- 
detic applications (see section 8). Precise measurement of the long 
term longitude accelerations on three operating synchronous satellites 
has already given an excellent discrimination of H22 H33 [8], as 

will be described shortly. The implementation of these effects to 
determine longitude variations in the earth *s gravity field is 
extremely simple. Aside from a rough determination of the initial in- 
clination of the orbit^ all that is necessary is to obtain good data 
on the geographic longitudes of the equator crossings of the near 
resonant satellite. In general^ the additional disturbances of the sun, 
moon and earth zonal gravity cannot be ignored where the longitude 
gravity effects (which are strongly dependent on inclination) are likely 
to be very weak. 
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DETERMINATION OF H22 AND H33 FROM 2 k HOUR SATELLITE DATA 

In the case of the 2 k hour satellites, numerical integration has 
shown that additional gravity perturbations to equation (67) are within 
the standard deviations of the nodal acceleration observations presented 
in Table iTS], This data has been derived from almost three years of 
tracking record on the free gravity drift of syncom 2 , J and the early 
bird satellite, all with nearly circular orbits. The gravity harmonics 
Table 1 are defined from the Xn,wx of ( 20 ) as: 

In ^ J 

SiVvNA Xiana^ 

•to conform to the common Way of expressing the series of longitude 
harmonics. The normalized longitude coefficients defined 

— — ^ ^ Tcrt-wV. 

L J j 

from which it can be shown that if all S are 1. the integral 

of the square of all the normalized harmonic functions over the unit 
sphere is 4 TT . The harmonics H22 an^ H33 reported in Table 1 were 

derived by a weighted least squares solution of (67) [more specifically^ 

\ 

using the orbit averaged coefficients of (33) in (l 9 ) withr = »i=/ ] ac- 
cording to the data in the table. The data was also tested for 
H42 and H44. These tests did not produce a significant change in the 
results for H02 and H33 or an improvement in their standard errors. It 
is concluded that the 24 hour data is not yet accurate enough or suf- 
ficiently widespread in longitude to allow a good determination of these 
other resonant harmonics through 4th order. 
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SENSITIVITY AND DISCRIMINATION OF LONG TERM HARMONIC EFFECTS ON NON- 
SYNCHRONOUS RESONANT ORBITS 

The fact is that only a limited number of harmonics are in long 
term resonance on a given commensurate orbit [i.e._, Eq. (65c)]. The 
possibility thus arises of using the long term evolution of these 
orbits to sense and discriminate these particilLar harmonics only. 

On a general orbit^ and considering the periodic perturbations 
of all the satellite's elements^ an infinity of longitude harmonics are 
active. In any gravitational analysis of limited satellite tracking 
data^ a major problem is the discrimination of one harmonic effect 
from another. The use of resonant orbits in satellite geodesy helps 
to solve this problem by reducing the set of gravity harmonics vhich 
are theoretically active. In addition^ resonance phenomena^ charac- 
terized by accentuated and persistant effects^ act as natural ampli- 
fiers of inherently very weak harmonic forces. Resonance increases 
many fold the sensitivity of these orbits to the weak perturbations, and 
thus enables the analyst to discriminate readily these weak forces 
from a background of much stronger (but not resonant) perturbation 
forces and observation noise. To judge the required sensitivity of 
the long term effects, we can use the accurate tracking period of the 
equator crossings of "early bird". This has shown [8] that we may 
expect to achieve resolution of the geographic acceleration of a near 
constant ground track to the order of 0.01x10 ^ rad. /sidereal day^ 
in about 30 synodic days of nodal tracking on the lower-r.ltitude sub- 
synchronous resonant satellites. With complete reduction of the 
relatively strong sun-moon effects of bi-weekly, monthly and longer 
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periods, we should be able to reduce this resolution further. Table 2 
gives a list of commensurate orbits of less than 2k hours and up to 
three sidereal days resonance period to illustrate the sensitivity 
and harmonic discrimination problem with respect to the likely long term 
effects. It is noted that for each resonant orbit except the 12 hour, 
the second relevant harmonic [from ( 65 c)] is two orders of n higher than 
the first. Since the strength of .each potential term falls off as , 

we can expect in general that the leading (lowest n) harmonic will 
also be dominant for that resonant orbit. However, the long term ef- 
fects are strongly dependent on inclination, with a different dependence 
for each harmonic. Thus it is necessary to consider a range of inclina- 
tions to determine both the sensitivity and discriminating powers 

of the resonant orbit with respect to the relevant harmonics. The 
evaluation in Table 2 is intended to be indicative of the likely order 
of magnitude of the resonant accelerations. It is based on the 
curve in Figure 8 which, in turn, is based on recent satellite-gravity 
results. The form of this "planning curve" was suggested by W.M. Kaula [l4], 
and appears to be a fair fit for the relevant resonant harmonics H 

through Hi 5 ,i 5 as far as is known to date. The normalized 
harmonics are defined from the non-normalized by: 



- 43 - 



(M 

w 

PQ 

C 

H 


I ° 

g ej to 

I 2 S 
g a ^ * 


s 'S 
2 1 •?: 
i&a 


^ O Q 

O cfi 
5 £, II > 


o R 

"S 


Q (D Cfl 




® R 
.^- ? 


Ic D 

K S W 


. « Q 


I -2 ^ 

I § I i 

■ m 6 B 
OS «i 


"tt .2 — 






^ooeo'-icv]cocoi30'4<«Ci-i«coot~t-''i*®i>'^*S£? 

«0C~N?0OO«5<N(NrtW0>e000^Tj<Q0«0£3^O 

rtt-l?OC<IOQO<Olrt^fOC<I-H.-IOOOU200^'TC^ 

1 eo N c 4 N i-< i-H iH rH .H t-H fH 1 -H .H rH in W « <N Irt 


) Tf< m <x! 00 03 o 1 


I m «D CO in t- < 




OOOOOOOOOO 

oooooooooo 


Mc5<oii5f^in^®ooooooo 
— — — -nt-ii-imoooooooo 
DOOOOOOOOOOO 




in o <M 
« in o 
en «£> t' 


<NCO|C')C-«e<NO'tt~ 
F-lrJ<O»O3C0Q0'^^*^ 

••t 00 eo 00 in t- 


ooloooooo 


in os 

Ol OS <o 

, N CO -It csi 

o in oj CO CD o 

‘ fh d d 


O (D I 


__~C-6-C~ 

in iH o N CO o 

o o o o o o 

o o o o o ^ 

o o o o o o 

till 


eqocDint-aoooMC'te'icoooosciicocviNFHOOBi 

c^TjiitF-iF^eOTjiincococoinoioosoooooO' 

05 OOf-iOC^Oi 1 *N 0 S 00 C 0 OCDtJ<OOOOOO( 

OOOOOOOOOOOC'JCOCDfHOOOOOOO 

dddddodooooooooooooooo 


I 


I 


-V 's. ^ 




i/ 5 oaoooMOC<ioiHinc<ic~c-oc-ooo 3 co»n 

F-isDca«oinc-inoo 50 scoiJ<cg«OTjiocoF-(OB 3 ?' 

inooO'ift>F-it~ososos«OTi<Tjico^oooooo’ 

F-IO'H(NOFHOOOOF-jFHCOC'JOqOOOOOO 

dF^ddddooooooooooooooo 


•-i -3- 'N. 

?o<oini>coooccjcMi-ioocoosTf<Ti<cocqiHOOMO 

t-eOC-COincO-<4*COCOC<JCCtF-iNFHOOOOOOOO 

-3<eOOFHOOOOOOOOOOOOOOOOOO 

oooooooooooooooooooooo 

ddddddoooooooooooooooi” 


inSeoF^ddoooooooocooooo olo 
t }«000000000000000000 o o 
ddooooooooooooooooo olo 


'^'*’'* 2 K 3353 W"^’* 1 ** 1 * • • • JWWS J 

K*’sTar 5 s"' “ - - 2 s 2 *2 2 ’ 

M ws go 2 2 2 if if -* 2 2 (5 S 2 2 * ^ S 

X X X S X X X X X XXXXXXXXXXX'. 


xxxxxxxxxxxxxxxxxxxxxx 


XKffiwwDsasasaswKasKwasacac: 


OOOOOCOOC-OOOO 

OOOOOOitO«STjiFHO 

doocOTji'dcoeoecieviC'jiM 


^ o o 

iHOOOOCOeOOOi 

t-oinocooocoo'^i 

F-lFniHcDOscDincso-^i 


( cq M CSI < 


'S .a 


© ^ 

X bc 

a § 

3 “ 

ii 

■a s 

O TO 


c o j 3 S 3 

e TO 

§ S. ® ■g 

•g « I -g 

TO *0 TO © 
«) O TO JS 

« X < U 

£86 


- 1^4 - 







As stated previously it can be shown [l6] that the normalized harmonics 
all have a mean square amplitude of unity over the surface of the geoid 
(actually a sphere of unit radius) if . The underlined values 

of Table 2 indicate those resonant orbits where both sensi- 

tivity (to the O.OlXlO"^ rad./sid.day^ levei in JO synodic days) and 
discrimination of a sihj^le harmonic effect appear to be good. Dis- 
crimination is considered good if the dominant (underlined) effect is 
likely to be an order of magnitude greater than the subdominant. It 
can be seen from Table 2 that only for the 8^ 12 and 16 hour satellites^ 
will 50 ^ inclined (or maximum payload) orbits be sufficiently sensitive 
and discriminating to harmonic effects in the above sense. (the ih.k 
hour orbit may be marginal) . Admittedly^ these tvo criteria are some- 
what arbitrary. One would like to use the already strongly discriminating 
properties of the resonant orbits to the maximum extent possible to be 
able to define unambiguously as many relevant harmonics as possible^ 
with a minimum number of observations. It is noted that for all the 
orbits in Table 2 except the 12 hour one an infinite series of sub- 
dominant harmonics have the same ^ or longitude frequency as the dominant 
harmonic. Thus^ investigating only a single resonant orbit (of constant 
inclination and major axis) will never allow an unambiguous discrimina- 
tion of these harmonics no matter how wide a longitude survey of this 
orbit is made (by natural libration or onboard propulsion longitude- 
shift manuevers). The technical reason is that there is no longitude 
discrimination between such constant harmonics in equation ( 67 ) vrhich 
governs the evolution of such orbits. 
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It has been proposed [17] that the gravitationally sensitive and 
discriminating properties of resonant orbits discussed here, be used to 
define specific longitude variations in the earth’s gravity field. If 
only one satellite of medium (30°) inclination is available for such 
purposes. Table 2 shows that H33, H32 and H43 may be detectable unambi- 
guously, in a reasonably short time, by simple single nodal longitude 
observations on 8, 12 and I6 hour satellites respectively. It is a 
property of the latitude dependenteof the potential function [in 
particular, in (65a)] that for the (dominant) 

resonances, the leading subdominant effect jn^ can always be "tuned out" 

(made zero) at one or more orbit inclinations above 30° ^ 31 } M- 
Such higher inclination resonant orbits may be particularly useful, 
geodetically. In addition, the inclination functions 
^n w)-fn increase in sensitivity to t as n increases. The effect 
is that "tuning out" of the subdominant to the dominant tesseral ef- 
fect ( is also possible at specific inclinations closer to 

30° than the first "tuning out" inclination for the dominant tesseral. 

Thus, for examplej a single geodetic 60° inclined satellite of 4, 

4 . 8 , 6 and 8 hour periods probably could unambiguously define H93, Hys, 

^55j H54 and H33 respectively in a reasonable length of time. It ap- 
pears from Table 2 [and a study of (65a)] that the strong resonances of 
higher order n (on orbits between I.5 and 3 hours) demand carefully tuned 
inclinations of greater than 60° to allow for sufficient one-orbit dis- 
crimination of effects. Of course, if other data on the harmonics can 
be used reliably as supplementary information, non (or weekly) discrimi- 
nating single orbit information may still be valuable. Thus it appears 
that all the one day resonances with periods greater than 1,6 hours 
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* • I? 

(reasonably drag free)^ may be geodetically useful if , We 

conclude with the observation from Table 2 that no resonant orbit of 
2 or more day synodic period seeds to have geodetically useful long term 
effects except the 16 and l4.40 hour satellites. These resonances are 
generally weaker than the one day ones, since they generally involve 
longer period, higher altitude orbits of high order effects. 
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LIBRATORY CHARACTERISTICS OF NEAR RESONANT ORBIT SATELLITES 


It should be clear from (67) that if the inclination and semi- 
major axis undergo only small changes, the long term longitude drift 
of a nearly constant ground track orbit due to any relevant harmonic 
^nm analogous to that of a circular pendulum. In fact many 

integrations of the complete eq.uations of motion have confirmed this 
assumption for a wide variety of resonant and near resonant orbits 
[9], fll], [12], As such, the drift as a function of time (more 
precisely the time as a function of drfit) can only be given exactly 
(and only for a single harmonic effect) in terms of an elliptic 
integral [9], Fl 2 ]. However, we can easily find the drift rate 
regime and also determine the critical or equilibrium points of the 
pendulum drift directly from (67)* Additional information about the 
libratory drift, easily found, will be the minimum oscillation period 
(about the stable equilibrium points) and the maximum drift rate pos- 
sible for "capture” in the libratory regime due to any relevant hamionic. 

We will take our standard pendulimi equation for the drift due to 

H in the form; 
nm 

where, for convenience, X refers to the ascending equator crossing 
longitude. Xg is seen to be a position of stable equilibrium ^ Ae^s') 
if ^ unstable equilibrium C Xe^u.) 



From ( 68 ) these positions are separated by J — Xe^u. ” * 

Equation (67) can be put in the form of (68) by rewriting (67) 






where 


6L / 


Now, since IT <1 0 , the positions of unstable equilibrium are given 
from (70) and (67) by writing: 

or, equivalently [using (63) with/=Z and the fact that n-m is an 
integer or zero]: 


Xe u. = X,„ 


9 

Similarly, those of stable equilibrium for given from 

(72) and (69) ty 

\ ^ JT -h XrtM ' 

, 0 « 

It may be noted that for low inclination satellites ( 
is always positive [3]^ and for these orbits (j2) and ( 75 ) te the 

equilibrium criteria. For example, for the 2k hour satellite the 
relevant all those for which n-m is even [see (65)] or 

:=: ZKj . Thus, for the low inclination 2k hour satel- 

lite, the unstable equilibriimi ascending equator crossing longitudes 
are at: 
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Xfe w- XiA»v\ "" 316 


ilh) 


and the stable equilibrium longitudes for low inclination 2 k hour 
satellites are at: 


Xe,5 ■" ArviM 


y 


(75) 


where C 


K- = = o \ ‘ 

CL / > / 


The situation is just the reverse for inclinations (at medium and 

I 

high inclinations for example) where these cases the 

stabl e equilibrium longitudes are given by (72) and the unstable by (75)* 
Thus : 


A 




X dy a- 


X*^ / 


= x„„ 


J 


(75a) 


(76) 


(77) 


/or resonant orbits where F <0 

• ' AM P 

Of course^ since the unstable and stable positions are equally spaced 
around the globe and separated by ^ radians_, there will be a total 
of of these around the equator. There will thus be m unstable 
longitudes and m stable longitudes in all for each relevant H 

YirA • 

The maximum libratory excursion of the ground track will thus be 560^ /m. 
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Proceeding with the integration of the pendulum equation (68), 
the variables tX) and \ can be separated by means of the dif- 
ferential reduction i 


\ = 


= iLxf = itxf 

iTTt a J.X 


First we rewrite (68) to apply to Xa = 5 ; 

• « X 

X I 5(aJ w(X-Xe^s') =■ ^ • 


(78) in (79) permits us to separate variables: 


JLC^] — — 2. I ^ X 


The integral of (80) is the energy or first, integral of the pendulum 
equation: 

( X ^ ^ I — Xe,s) . 


If the initial condition X = X^^at X=x« is given^ then (8l) may 
be written as: 
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From (82O we can easily find the maximum drift rate change possible in 
a resonant orbit due to each relevant harmonic. This will occur with 
nearly exact resonance ( X^-sso) momentarily over a position of in- 
stable equilibrium. Then the drift rate passing the stable point 
radians away from Ae,t4. ^ will be the maximum possible 
libratory drift. 

For ~ ~^/i^ , (82) reduces to: 

The maximum libratory drift rate (a.t X=Xe^s) will then be : 

( \) =z ^ , 


Since the maximum acceleration in the libratory regime, from (79) 


IS 


j ^ ^h.h1 • we can write (83) as: 


(-X) 


fA A Ay n ivv 


= z. 




hn 


Another parameter of interest with respect to the regime of the 
resonant satellite is the minimum libration period. This is the 
period of small "pendulum" oscillations of the ground track about the 
stable equilibrium points defined in (j 3 ) ^7^)- is 

always sufficiently small (as in a simple pendulum oscillation)^ 

can be replaced by and (79) becomes a 


simple harmonic equation 


From ( 85 ) we can easily find the small oscillations ot minimum period 
libration about the stable equilibrium point. This minimum period is: 


-r 


Mm 


Xir 

LAj 


^ a.ir 
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Table 3 below gives these 3 libratory characteristics 
and for the most geodetically promising (underlined) orbits in 

Table 2. The X . coefficients in this evaluation are taken from the 
"average ^ curve of Figure 8 and not from any specific set of coefficients. 

A comparative case, from actual data on 2k hour satellites, is provided 
at the bottom of Table 3 . 

As mensioned previously, the complete analytic treatment of the 
drift regime due to a single resonant harmonic (for a circular orbit) 
can be carried through by solving (82) in terms of elliptic integrals. 

This will not be done here. The form of the general solution has been 
shown previously the specific solution for the long period 

of libration of one day resonant orbits [ 3 ], [l83. As in the circular 
pendulum, the actual libration period is strongly dependent on the libra- 
tional amplitude. In addition, Gideon et al. (I 966 ) [18] has also obtained 
complete analytic solutions for the one day librational resonances of an 
eccentric orbit due to single harmonics. 

We note only that since the basic averaged equation of librational 
motion, (6j)^ is nonlinear, the complete solution of the quasi-libration 
due to the total, infinite set of resonant harmonics, cannot be a simple 
superposition of the individual harmonic solutions . 


CHARACTERISTICS OF GEODETICALLY PROMISING RESONANT CIRCULAR ORBITS 



♦From actual data 




In fact there is no second integral to (6'j) which is known in a 

simple closed form. However the first or energy integral to (6'j) 
can be simply written down analogously to (82). In terms of or 

the longitudes of the ascending equator crossing, (6j) may be rewritten 
with the differential reduction ( 78 )^ Q-s: 

xcx) - 2.ixn ^ C.OS S i^j 

relev/c^-n*t 



The integral of ( 87 ) in terms of the initial conditions 

X = X(,A.t \ = X^C±-c>) is: 


Then, formally, we can separate variables in (88) and find the time as 
a function of drift from: 

X 

— j 1 X j ( xn fsi/^ (m[\- Xrt,v)^ -f- S ^ , / 

' J i re\. '' ^ 


Of course^in all of this simple derivation, we have assumed the elements 
of the orbit to be essentially constant (in particular only 

over the resonance period n^ but over the long librational time t. 

While this will undoubtedly be a very good ass-umption for the semimajor 
axis (which in most cases will suffer small relative changes over long 
periods of time), it may not be so with regards to the inclination. 

Due to the gravitational action of the sun and moon combined with earth 
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zonal gravity^ the orhit plane of a satellite (otherwise imperturbed) will, 
precess about a plane between the ecliptic and the equator causing a long 
period sinusoidal like change in the satellites inclination [isl , The 
period of this inclination change ( with an amplitude of the order of 
about 10°) can be as low as 10 years for close earth satellites. However, 
for the close earth satellites ( ki 4~ , for example) zonal gravity 

equatorial precession predominates so that the total inclination changt; 
is small. Thus, it will probably be found that only the relatively 
weak far earth resonances can be seriously disturbed by these plane 
precession effects. 
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DISCUSSION 


First, the use of the term "resonance'' and "resonant orbit" in 
this report should be made clear. In mechanical systems resonant 
phenomena occur when the periodicity of applied loads are close to the 
natural (unforced) periodicities of the system. Resonance in this 
sense is characterized by greatly amplified vibrations of the system. 
Viewed internally, the applied loads are always in step, or in the same 
direction as the internal parts upon which they act and so continvially 
increase the vibrational energy of those parts. But viewed externally, 
the loads are purely sinusoidal and so the mean translational energy of 
the center of mass of the system is unchanged. In the case of the 
satellite whose orbit yields a constant ground track, we have seen 
that, with respect to certain harmonics, the applied loads are not 
purely sinusoidal but contain a bias component which drives the long 
period libration. In the case of the constant groijind track orbit, 
what we have called resonance implies a change in the translational 
or orbit energy of the satellite . The vibrational, or internal energy 
of the orbit can perhaps best be characterized by the eccentricity. 

In an eccentric orbit the satellite can be imagined as vibrating con- 
tinually about its mean anomaly and semimajor axis, or mean radius. 

In this report we have chosen to ignore this vibrational aspect of 
the orbit by orbit averaging the disturbing force from the outset. 

We do this because we are here only interested in long term orbit 
energy changes. In terms of the mechanical resonance analogy, the 
full effect would be equivalent to a mass-spring assembly on wheels 
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vibrating and translating under the action of a biased, sinusoidally 
varying external force. We have only concerned ourselves with the 
steady, translational, part of the complete motion. In celestial 
mechanics these changes are called seciilar, or effects which increase in time 
without apparent limit. In the case of commensurate orbits, the bias 
of the disturbing force will also change sinusoidally over a long 
period of time. This is a phenomenon which we call libration, having 
as its proper analogy the circular pendulum. 

Resonant orbits as they are used here are synonomous with 
commensurate or constant gro\md track orbits; those whose rates are 
rational with respect to the earth *s rotation rate. In addition, of 
course, we have specified the orbits to be near circular. The '^resonance 
phenomenon that we have, by our orbit averaging method, limited to our- 
selves to, might better be called secular or librational-circular pen- 
dulum effects. (In particular we have only really dealt with the secu- 
lar changes in the semimajor axis). But the "in step action" of the 
longitude harmonics on these commensurate orbits give rise also to 
true vibrational resonances ( in radius and anomaly) following the 
mechanical analogy [ij]^ [20], [21]. Therefore, though the effect we 
analyse for should not commonly be thought of as a resonance effect, it 
does arise most strongly for those orbits which do show true "resonance" 
phenomenon and thus these orbits may be justifiably called resonant 
orbits. But when we speak of a "resonant gravity harmonic" in this 
paper, we use the worii resonant only in the sense of being capable of 
producing the amplified changes of libration which are also properties 
of the resonant orbits. 
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It may be pointed out that the true vibrational resonances in a com- 
mensurate orbit imply a buildup^ or amplification, of the orbit ec- 
centricity. This amplification would be without bovind if true commen- 
surability could be maintained. Of course it cannot, since the resonant 
orbit generally librates. In fact, as might be inferred from its depen- 
dence on the direction of the dist\irbing force, the eccentricity itself 
in a resonant orbit appears to go throu^ the same long period libration 
that the semimajor axis does [l8]. However, these true "resonances of 
eccentricity" may not be as useful geodetically since all the longitude 
harmonics H appear to contribute to them. 

Finally, we would point out that the pendulum like changes of the 
semimajor axis [describable from a seaprated solution of the coupled 
equations (66) and (67) ^ see Ref. [9] for example] have a periodicity of 
the order of years only for near commensurate conditions . These are the 
librating orbits first discussed in Ref. 10 for the 2 k hour case. In 
cases where the ground track moves at a considerable rate, the basic 
circular pendulum equations [(66) and (67)] with small modification 
(along the lines outlined in Refs. 3 and I8) may still serve as the 
model of the evolution. In these cases of fast global circulation of 
the ground track, evaluation of (66) and (67) shows that the crossing 
rate and semimajor axis oscilate with a more rapid frequency but much 
diminished amplitude as compared to a libration case. For example, in a 
gast circulation case the period is typically of the order of days and 
the amplitude may be of the order of hundreds of meters ( for the semimajor 
axis change) . In full libration cases well away from the equilibrium posi- 
tions, the amplitudes of the semimajor axis change may be of the order of 
tens of kilometers T9]. The fast world circulation regime (far from resonance) 
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has been called dynamic resonance by Blitzer [2l]. But by the dis- 
cussion here it is actually a circulating pendulum and not a true 
resonance phenomenon. The reduction in amplitude of the pendulum 
effect in a fast circulation regime is due to the rapidity with which 
the orbit averaged bias force averages out over a global circulation. 

But in spite of the small amplitudes, "circulating penduliim resonance" 
has been used to good effect in discriminating a number of very high order 
coefficients in the earth's field from a dense global tracking of the 
Navy's Doppler-Transity Satellites [I5]. 

The limitations of the model used here to derive the key orbit evolu- 
tion equations (66) and (67) should be re-emphasized. Though strictly 
speaking these equations apply only to circixLar orbits of exact resonant 
period, they have been shown, by mmierous examples [8] and in a rigorous 
■way [3], [18], to be applicable with only slight correction or modifi- 

cation to a much wider class of satellites. This class includes those 
with drift rates considerably exceeding the maximum permissible for 
libration as well as for noncircularorbits with eccentricities no 
higher than about 0.01. The equations, furthermore, give only the gross 
orbit averaged effects due to the earth's longitude gravity field. But 
over a period of the order of month's, it may be presumed that, well 
away from any strong influence of the earth's atmosphere or the moon's 
gravity, these equations will describe the dominating long term effect 
on the geographic configuration of all near resonant near circular 
orbits with the exception of Inclination change due to the sun and moon. 

The extension of these regime equations to the class of eccentric 
resonant orbits is indicated by Gedeon, et. al. [l8]. 
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SUMMARY AM) CONCLUSIONS 


The long term geographic evolution of circular orbit satellite 
ground tracks which are originally near stationary, has been found. 

The eq.uation governing the motion of the near stationary groimd 
track is essentially that of a circular pendtilum for each relevant 
earth longitude harmonic. The relevant harmonics for circular resonant 
orbits form a tenuous but infinite subset of the infinite set of longi- 
tude harmonics 

For orbit periods very close to resonant, or rational with respect 
to the earth's rotation rate, the sinusoidal ground track librates with 
a maximum excursion of 560°/m due to an individual relevant 
(resonant) earth longitude harmonic. The complete period of libration 
is the order of 2-10 years for the strongest resonances and depends 
strongly on the libration amplitude and inclination. 

These libratory-like evolution equations can be extended to far- 
from-stationary conditions to cover the smaller^ but detectable^ drift 
oscillations of higher frequency, called "dynamic resonance" by Blitzer [2l]. 

Circular resonant orbits of periods 3, k, 4.8, 6, 8, 12, l4.4, 
and 16 ho\irs appear to be particularly suited to discriminating 
mambiguously the longitude harmonics Hgs, Hyq, H55, H54, H33, H32, 

H55 and H43 respectively, in a reasonably short period of time. The 24 
hour satellites have already provided the best discrimination of the two 
leading sectorial harmonics H22 and H33 as of I966. It appears that a 
single controllable test satellite, "hopping" from low to high circular 
resonant orbits and "sliding" a number of times within each orbit could 
provide unique discrimination of most of the above harmonics in a 
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year or two. The orbit determination for the resonant satellite used 


as a gravity probe would not have to be elaborate . It would have to 
insure that a reasonably circular orbit has been achieved. The inclina- 
tion must also be determined to about 0,1° accuracy. During free drift 
periods the longitudes of all the equator crossings is all the data 
that is necessary for a rapid harmonic determination. 
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APPENDIX A - THE RESONANCE SELECTOR FOR SUPERSINCHRONOUS ORBITS 
Introduction 

It may be verified that, for the supersynchronous k 8 hour satellite 
[since Y~ ~Tj_ =■ My and w = y~ s i for the i<-8 ho\rr 

satellite] the half orbit periodicity conditions (29a, b, c) suggest 
that the resonant harmonics are. 

For n-m even: 

m = 2, 4 , 6, . . . 

For n-m odd: 

® 3, 5^ . . . , or 

^ 42 ^ ^62 • • • y H44, H04, , ,, ^ 4 -Z.y , ,, 

^43^ Hea, Hs3, . . . 

However, criteria i'^ 2 ), extended from Allen’s work [ 3 ] and apparently 
equivalent to (29a, b, c) for eliminates some of the leading 

harmonics in the above series . Thus, (32) for the 48 hour satellite 
is: 

2m( relevant) = n( relevant) — /a — '/ y 

which produces the following table of relevant harmonics: 

m n 

1 4 , 6, 8, ... 

2 4 , 6, 8, ... 

3 6 , 8 , 10 , ... 

4 8, 10, 12, ... 

For example, criteria (32) says that H22 is not resonant on the 48 hour 
satellite white criteria (29a, b, c) says that it is. 


- 64 - 


We will now attempt to verify by direct calculation, which of these 
two resonance criteria is correct for the hour satellite and pre- 
sumably for all supers ynchronous orbits. Equation (B8) in (B6) gives 
the along track force on the 48 hour satellite C h ^') < 3 .ue 

to H22 as: 


F 

=r S|a/2-(Ao-X 

k; 


2.2. 


\cos&^ 

[ 5 ,-^ 8 , cco^i 

2 ^ *^L 1 


CCS^?< 


+ [, - > w-g, ( I 


r il r 




57 ^ 


lc£.#©j 



r J 

J ^ 



i' 

(Al) 


Note that f in (A|) has been written in terms of for the 

48 hour satellite. The orbit average of (A|.) is to be taken over 
ZTT since^ for the 4-8 hour satellite, r = 1, 

It can be seen from (Al) that the terms in S 1/^ 2- ^ Vo are all 

odd functions with respect to ^5 ^^ 1 % ‘ Similarly, the tenns 

in C05 T^LXo-X'vt?^ all odd functions with respect to £% =TT 

Thus the orbit average of for a 4-8 hour satellite is 

zero and H22 is not resonant on this orbit. Thus, even though 

^^/■2-i^^4-8 hour) is periodic over a half orbit [ from (29b)] it 
orbit averages to zero. This can only be the case if it also half-orbit 
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averages to zero. It is interesting, but so far unexplained, that no 
half orbit periodicity in Fy has been found to orbit average to zero 
for subsynchronous resonant orbits. For resonant orbits higher than 
synchronous, however, we presume from the inference of the preceeding 
calciilation that many such periodicities will half-orbit average to 
zero (in the manner of products of periodic orthogonal functions). In 
particular, the inference is strong that (32) is still the correct 
resonant harmonic selector for supersynchronous orbits. 
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APPENDIX B - ALONG TRACK FORCES FROM LONGITUDE GRAVITY HARMONICS THROUGH 
4th order, for RESONANT ORBITS 

It will be useful, in discussing the orbit averaged behavior of 
a resonant orbit vinder the action of a longitude gravity harmonic, to 
develope the along track force from this harmonic as a fxmction o'f 
earth angle turned since the satellite was at its ascending node. 

Figure 2 illustrates the geometry applicable to the consideration here 
of a circular orbit satellite. From this figure we note the following 
spherical trigonometric relationships: 


5(^/<X COSi-^coS^ 

CO 5 o<. — ^ /-TAti 

<p — 

TAiOt^s) 

The geographic longitude satellite at any time t after its nodal 
passage is: 


(Bl) 


I X — X o ^ L ” • 


(B2) 


We note from (2) that for the resonant orbit specified by whole posi- 
tive miners and > (B2) becomes : 


X = X. V }j^) - 1^1 + . 


(B3) 
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We also note that for the resonant orbit can be replaceyjin (Bl) 

by . The along track force on the circular orbit 

satellite, due to H is (See Figure 5 ) ; 

’ run o >/ • 


F- 




C.OJ5 o(^ 


(Bh) 


The longitude and latitude components of the earth gravity perturbation 
forces through 4th order are given in Section VIII of Ref, [ 6 l], It 
can be seen from (B1 ),(b 4) and Ref. [Bl] that ve can write 


F 


Tj 




(B5) 


where is a constant over the circular orbit. Substituting 

.i 

' COSL TA^jtfie ) 

f7 ^ 


(B3) into (B 5 ), with TAf^ ( CoSt-Xhtit^^e. ) gives: 
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The g's are all functions of i, n, m, and t}i^ only: 

M iX 

COSI^SX 
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For the 2 h hour satellite and the g’s can be evaluated 

through 4 th order (for Wf-o) from equations (31), (36), ( 42 ), ( 45 ), 
(54), (60) and (63) in Reference [B 2 ], On the right sides of these 
equations we replace 0 by and find (for Hxx) : 


^ ^ I- VO 




- Cos 




If,. 

it 




_ 5 


^ i 5VV 


6itJzC^-^)(^osu w) 


a := 




— 5 itJ z!^^) (\ - C,osi) 
i* ^ 




(b8) 


For H31, we find; 
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For H32, we find; 
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For H 33 we find; 
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For H 41 , we find; 
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For H42 we find; 
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For H43, we find; 
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» For H 44 , we find; 
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APPENDIK C - THE THEORY OF THE 12 HOUR SATELLITE 


We can derive the theory of the 12 hour (circular orbit) 
resonant satellite in two ways. Essentially they are both methods of 
evaluating the orbit averaged disturbing force. In the first we com- 
pute directly the average along track disturbing force from the associ- 
ated Legendre harmonics as they appear in spherical coordinates 
(Appendix B) . In the second, we compute the potential as a function 
of the Keplerian elements (see Section 5 )- Th® derivative of this with 
respect to the orbit argument also represents the along track force 
whose orbit average may be taken. Once the orbit averaged force is 
obtained. Equations (l8) and (19) give the orbit evolution as long as 
conditions are reasonably close to resonance. The first method was used 
in Reference [Cl] to derive the theory for the 2 h hour satellite through 
4 th order in earth gravity. The results of the orbit-averaged compu- 
tation appear in (jU). These may be compared with the results of the 
second method of computation which is summarized in equations (66) and 
(67) for any resonant orbit. The two methods give identical results (as 
they should) for the 2 k hour satellite 

It is noted that the straightforward method; , while cumbersome to 
carry out, is capable of generating the orbit averaged force for any har- 
monic-resonance condition through evaluation of the ^ equations in 
Appendix B to any order. However, extensive evaluation of these equa- 
tions for synchronous and nonsynchronous resonant orbits convinces us 
that calculation of the orbit averaged force through (65) will yield 
correct results in all resonance cases. In particular, where the reso- 
nant harmonic selector. Equation (52), has been violated for many har- 
monic-resonance conditions at, below^and above synchronous altitudes. 
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^ the straightforward orbit averaged force, as given throixgh i»-th order in 
Appendix B, has always been zero. 

We show here, for example, the correctness of (65) for evaluating 
the H 32 (12 ho\ir) resonant effect. In Appendix D we show by numerical 
integration of a partical trajectory that ( 67 ) is correct for a speci- 
fic configuration of a I6 ho\ir orbit Cr=*3^ W . This latter case 

is important because it represents a check of the extension here of the 
formulae [similar to (66) and (67)] given by Allan [3] for the n*=l 
day resonances. 

Before beginning the evaluation of the leading resonant harmonic' 

effect, we may note that in the case of the 12 hotir satellite, all 

the nonresonant harmonics appear to be selected by the "antiperiodic" 

criteria ( 50) and ( 51) . For example, the strictly nonresonant harmonics 

for the 12 hour orbit = from (jO) and (31)- 

H , where 
nm 

m = i<-, 8, 12, . . . for n-m odd, 
m = 2, 6, 10, . . . for n-m even 
and m = 1, 3, 5W 

It is seen^then^that from these strict negative criteria, the harmonics 

which may have long term resonance effects on the 12 hour satellite 

vare those H Where 
nm 

m = 4, 8, 12, . . . for n-m even, 
and m = 2, 6, 10, . . . for n-m odd. 

Indeed, such harmonics satisfy the periodicity selector (29) which, 
we have seen^ is strictly equivalent to the more general selector ( 32) 
or (65c)^ for l^>M*( subsynchronous resonant orbits). 
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Thus^ the lowe.st order (n) resonant harmonic on the 12 hour orbit 

-W 

is H32. It is presumed to be dominant both because of the poten- 

tial decline -with distance [See Equation (67)] and the apparent 
decline vith increasing n. The normalized coefficients themselves 
give a correct interpretation of the mean realtive strengths of the 
harmonics at the surface of the earth. It should be noted^ though^ 

(see Table 2 , for example) that the dominance of harmonic effects is 
strongly dependent on inclination as well. This dependence reflects 
the integrated effect of the satellite's motion in latitude. 

The spherical harmonics themselves^ of course^ vary strongly in 
latitude.. When we speak of gaieral dominance here^ we refer to the 
mean effect over- all inclinations . 

Proceeding with the evaluation of the long term effect due to H32 
(12 hour)_, from Appendix Equations (B6) and (BIO)^ we have (with 

a 9 from Eq. (h) for the 12 hour orbit); 

1 ^^ J 




i 2. f f - 

Cos I 

C I— 3 ; 


- 5/^^26't. Cecil -|)j 4 - 5 

sll -h C^Cci2i\c-^3i^ / Sl>^ 


J i- — — 


(Cl) 
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It can be seen that the factor of above is an "odd" 

fmction with respect to ©e — . Thus, since the double-orbit 

average of (Cl) [See Eq. ( 9 )] is achieved in the interval 
in — i sidereal day) , the term will double orbit average 

to zero. Retaining only the relevant CoS term of (Cl): 


2.0-5/A;'’;2^e^/A^X) 


(C2) 


Let ; 

A = s,fj\e^-i) -i- 5 /a< z©e ‘ 

Then (C2) becomes: 


-h Slhf\&e -/) I i -^/A» 2<?e (» )j’ - Co^ze^tufJu [ 3>/vjt5>V 2^>£-/)£ 


5(>^ 4^ } 






U-c^si)s,yi.s,^4-ee Lxs,^\&^ oosl( 
Hi- 5/AJ^-2.e^ ) L 


« Co^XB^ ^,^L-i)i \-cosi )~j ? ^ 


(C3) 
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Let B be the term in [ ] immediately above. 


Expanding B, ve have: 

^ i\-Lo$V) t za^Lt-cofi J 

- / 4- Cfisl) ^ZCOfC -t 2 . ceil 20^ 6 1 - Co^l ') 

~~5c.o^z&^ ^ co^ze^i\~c.os i ) » 


Red-ucing B above: 

B— — 5//J Z6^L I - coil. ) ^ Co^ZG^i, l~caSt. ^ 

^ I- Coi i) ( Co^ Z6g_-i- Z0^) — / —2.Z-C>SL^» 


Further reducing B: 


B = 4 6i-co5a')- l-±cosi 


= - 3 ^.o 5 4 C I - $ ©e 5>/J V ) . 

Combining B (above) and k, (C3) becomes: 

+ - i£f;£f- i,\-c.05i)ls,iii 5>/jV^e) f » 

^ (C4) 


In the 12 hour orbit, ^ [see Equation (la)]. Thus, 

a.Xe^ and Equation ( 9 ) becomes: 

(for the 12 hour resonant orbit, double orbit averaged): 


(^02. Hou/J.) 


a.*iT- 

J' ' 


(C5) 


We use the following integrals : 

air ° ° xyr 


2.1T 


i 

anr 


JL II >1. II 

I Cc^ J ^ 0e = 
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- i. ^;a> 4^ Co546e = ^ ft') 


4 


_ I 


g 


2.7T 


\ 

aTT 


ivtr 

^ir 


o 

XTT 


t 

3. 


- _L 1©^ H-5/^i2.£?eC^>5Z»el - J_ . 

J * ^ -4-TT ^ — 7 X 

JLTT 

[ 1 % - 
f = fS- 1 ^ . 

<J 

With these integrals^ (Cil-) in (C5) gives: 

r = |c, i-o>;2./Ae-X:!,) -5/^^^* ^^,tiUo6L 

\iX ~Q— 

-J_ Co5li\-Co$i jsw.;] 


■)l] 


t-COSi. 


^ [^ I - ZC0$u -Scos^cj . 


(C6) 

Equation (C6) in (19) gives the long term acceleration of the ascending equator 

i 

crossings of the 12 hour satellite due to H32. We note that F= F<^ 

Or 

and X— Xo unrestricted position of the crossing over many or- 

bits^ treated^ effectively^ as a continuous variable. Thus^ this formu- 
lation gives : 

I<1» 

(C7) 

We now evaluate Equation (65) fo^ "the second formulation of this ac- 
celeration. Evaluating (65c) for n(relevant) = 3? m(relevant) = 2 , 

Jl- (12 ho\ir) = 1/2: 


- 79 - 



p (H32, 12 hour resonance) — ^ ~ I 


Thus, (65a) gives the resonant inclination function for 1132(12 hour) as: 

* * |c 

^ 1^2. Co5^Ct/^)^/jJLy^) - ' 


But^ 


Thus: 


5jA>4i/o.) = ^ ^ c<5 scl/^) - ^i/^titcoslj 


f'li) =(f)[ 

55./ 


A 

- 4: i\irCo$i. )‘J- U-cosC 


^ JJl C\^ce>5l')^ 1^- ^05/ -X Cl-coi<. 

65 

== 'lElUlt ^ \ ^ xcosi -Scoi^z^ ^ 


From Equation (C8) in (67) ^ "the acceleration of the crossing longitude 
due to H32 in the 12 hour orbit is • 


X =r XA-ir' 


j - '.LLSih}i[\~XcoiC- SC.05 V) co^ 5. ^ X -X 3 a.^ S 

^ 6 ~J 


rad./sid.day^ 


which is identical to the direct orbit averaging formulation. Equation (C7). 
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Similarly, from either the direct orbit averagining formulation, or 
from Equation (67), we have found that the effect due to on the 12 
hour satellite is: 


X 


(due to H44 only) 


-7.4-ir 




rad. 

Sid. days 


(CIO) 


Tables D2 and D5 give a comparison of the theoretical motion due to 
H32 and H44 from equations (C9) and (ClO)^with numerically integrated 
partical trajectories. 
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APPENDIX D - ONE AND TWO DAY RESONANCE TESTS: 8, 12 AND l6 HOUR ORBITS 

We can check the validity of the formulation of the orbit averaged 
librational equations of motion (67) by comparing them with a complete 
numerical solution of the equations of motion for a trajectory in the 
presence of the critical disturbing forces. Tables D1-D6, show the 
results of these numerical studies on near circular resonant orbits of 
8, 12 and 16 hours. 


The trajectories run about I6 sidereal days from near stationary 

ground track conditions. We hope to determine that even in this short 

time span, it will be possible to discriminate the relevant harmonic 

accelerations to reasonably high accviracy according to the simple theory 

presented here. 

8 HOUR ORBIT 

i)'or the 8 hour (one day) resonant orbit the leading relevant 
harmonic is H33 (having the least n satisfying (52) for = ). 

The orbit averaged equation for the acceleration of the longitude of 
the ascending equator crossing due to H33 on an 8 hour orib [(67)^ with 
JL^Hj p ^ I ] gives: 

rad./sid.day^ ^ (D1) 


The normalized are given from the non-normalized [ in (67)] 


l^iM 




A fPk 


r in- 

i' C ] 


-m')! 




(D2) 
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Evaluating (Dl) for: "31 = | 4 - 7 f ^ ^ ^8 hour orbit) ss ^.fL. 

L the theoretical longitudinal 

acceleration in the neighborhood of \~0 should be: 


\ ^ 0.4'44'7aIo'^ ro.L!ulX^^. 

(D 5 ) 

( see Table Dl) . 

The following data is from a partical trajectory [generated by a 
modified Encke method [Dl]] computed numerically^ in the presence of 
perturbations due to H33 [with constants given above] and H31 only. 

The specifications for H22 in the trajectory generator are: 

Xxx specifications for H31 are: 3 ^. :rr X3, — 

From the theory^ H31 and H22 should have no long term effect on the 8 
hour satellite. By long term here^ we mean over periods which are 
multiples of one sidereal day. 


Time 

(Hrs.)^ 

T 

Ascending Equator 
Crossing Longitude 
(Degrees East of 
Greenwich)^ X 

Semimajor 

Axis 

(Earth Radii) 

Inclination 

(Degrees) 

Eccen- ; 
tricity 

(lO“s) 

0 . 

54.0 

3.1781 

30.0 

o.i 4 

191.476 

54.00290 

3.1781 

30.0 

0.13 

382.950 

54.02279 

3.1781 

30.0 

0.15 


TABLE Dl 

Numerical Trajectory Data for an 8 Hour Earth Satellite Distiorbed by '^31 

and H33 Only. 


In the above trajectory^ the longitude acceleration is given approximately 


as: 


» ■ 

X ^ 

M£A} 


I 

( 4 T)^ 




(d4) 
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where 


constant 


AT=r -T; — 

j ri 


Since ^ -p ^ 8 sidereal days, the "measirred" accelera 

tion in the numerical trajectory is: 


t: =!- 


rad 


STXU 


s±d..X 




The discrepancy between the orbit averaged theory and the numerical 
trajectory data is only 0.0014 XlO ^ rad./sid.day^ out of 0.4633x10 ^ 
rad/sid.day^ or about 0.3^» 


12 HOUR ORBIT 


For the 12 hour (one day) resonant orbit, the leading relevant 
harmonic is H 32 (having the least n satisfying (32) for f- is. y\*s i 
The orbit averaged eq.uation for the acceleration of the longitude of 
the ascending eq.uator crossing due to H 32 on a 12 hour orbit [(67) ^ 
with/sA p—l ] gives: 

^5^2. Ll-XCosl -3co^i) Co<,xL\~^3:^ 

rad./sid,day^ • 


X - 


__ iS'l.joX 


< 


y 


) 


We consider a partical trajectory with only H 32 acting as a pertiurbation. 
The initial circular orbit has a period of 12 hours ( = 4.16449 E.R.) 

The inclination is 30 ° and the initial longitude of the ascending equator 
crossing is at 0°. The Has constants are: — — \0 ^ ~ ^ 

From (D6) the long term acceleration should be: 

X _ ^0.34-61x10 

TH. rad./sid.day^ 


(see Table D2) 


■ The following ascending equator crossing data is from another numerically 
integrated particle trajectory, in the presence of perturbations due to 
H32 only (with constants given above). 


Time 

(Hrs.) 

Longitude 

(Degs.) 

(E.R.) 

i 

(Deg's) 

Eccentricity: 

0. 

0. 

14-. 16^4-5 

30.0 

10"® 

I9i.i4.76 

-0.00660 

i4-.l645 

30.0 

CO 

0 
1 — 1 

382.953 

-0.02598 

4.1645 

30.0 

10"® 


TABLE D2 


Numerical Trajectory Data for A 12 Hour Earth Satellite Disturbed By 
H32 Only. 


In the above trajectory the longitude acceleration (over the 16 sidereal 


day period) is given approximately by (d 4) or: 


\ =-L 1-0.02.538 +• 0 . 0 / 320 I /, = _o.34.B5xjo'-^ 

rad./sid.day^ 


Once again the discrepancy between the orbit averaged theory and the 


numerical trajectory data is very smal].^ amounting to only about 0.1^. 


The simulated trajectory summarized in Table DJ has initial condi- 


tions precisely the same as in Table D2 except for the addition of earth 


zonal gravity and the attractions of the sun^ moon and planets. The 
zonal gravity constants used are J20 = 1.0823X10 J30 = -2.3x10 ^ 

and J40 = -1.8x10“^. 


TIME 

LONG. 


i 


( Hrs . from I966 . 0 ) 

(Degs.) 

(E.R.) 

(Degs.) 

e 

0. 

0. 

4.1645 

30.00° 

10“^ 

191.434 

0.14563 

4.1645 

30.02 

1 

0 

t — 1 

382.870 

0.27110 

4.1645 

30.02 

io“^ 


TABLE D3 - Numerical Trajectory Data for A 12 Hour Earth Satellite 
Disturbed by H32, Earth Zonal Gravity^ the sun^ moon and Planets 
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In the above trajectory, the longitude acceleration is given approxi- 
mately as: 

x' 

64 - 

rad. /sid.day^. 


^D.llUO 


— -i).5'4-3xlo 




(D8) 


It is apparent from the 8 hour study that contributions to error in the 
simple theory from non resonant earth longitude gravity are negligible 
for these nearly stationary orbits. But from the previous two tables 
the error in the long term theory due to other gravitational effects 
is clearly not negligible over 16 sidereal days. In fact, other studies 
have shown^[D2], [D5]'|that the sun and moon alone accounts for nearly 
all of this major discrepancy. But the presence of the strong nonresonant gravi- 
tational effects alone gives the particle trajectory a mean motion over 
only 16 days close to that of the fully perturbed motion. Therefore, 
it appears reasonable to try and evaluate the "sun and moon" effect 
independently of the fixLl trajectory. Then, subtracting this from the 
full effect, we hope to get a reduced acceleration which can be assumed 
to be due to the relevant resonant gravity terms alone. 

The simulated trajectories in Table are a test of this error 
reduction idea. Initial conditions were the same as in Table D3 
except ^ earth longitude gravity effects were included for the tra- 
jectory with bracketed data. The rimning time for both fixLl and non- 
resonant trajectories was l4 days. 
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time 

LONGITUDE 


i 

e 

(Hrs. from I966.O) 

(Degrees) 

(E.R.) 



0. 

0. 

4.1645 

30.0° 

10"^ 

(0.) 

(0.) 

( 4 . 1645 ) 

(30.0°) 

(lO"^) 

167.505 

0.12756 

4.1645 

50.0° 

10"^ 

(167.505) 

(0.13243) 

( 4 . 1645 ) 

(30.0°) 

H 

0 

1 

355.011 

0.24023 

4.1645 

30.0° 

lO"^ 

(355.010) 

(0.25974) 

( 4 . 1645 ) 

(30.0°) 

(io“^) 


TABLE Numerical Trajectory Data in Two Sun and Moon Perturbed 12 

Hour Orbits; Disturbed (unbracketed data) and Undisturbed 
(bracketed data) by H32. 


In the H32-disturbed trajectory (unbracketed data)^ the measured ac- 
celeration is: 


\ 




(fully disturbed) 


4-9 ^ 










-r 


rad./sid.day^ 


(D 9 ) 


In the H32-undisturbed trajectory (bracketed data), the measured ac- 
celeration is: 


X (partially disturbed) 

M£A5. 


^_L [o.XS374--O.Z(,^el]/sy.2.3i, 
4-9 / 


= — rad./sid.da^'^ « 


(DlO) 
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Equation (DlO) subtracted from (D9) gives the reduced^ or sun and moon 
corrected acceleration as: 

X (corrected) = -0.3^76x10 ^ rad./sid.day < (Dll) 

But the theoretical acceleration due to H32 alone for this orbit is 
very close to (D6). Therefore ve see that this simple method for evalua- 
ting the earth resonance acceleration is valid to within about 0.15^ 
with this particular orbit. . JFor relatively short trajectories^ where 
the long period mean motion is still approximated very well by the 
nonresonant gravity effects^ this simple correction technique should 
apply quite well in analyzing resonant orbit data. For longer tra- 
jectories^ it may be necessary to alter the initial energy of the simu- 
lated trajectory without resonan<t gravity so it approximates the fully per- 
turbed trajectory more closely. This may be done most simply by altering 
the initial semimajor axis^ for example. At any rate, tests of the correction 
method, similar to that above, shoirld always be performed to check the 
method. We may note that this is, if successful, a direct method for 
evaluating the long term contribution of nonresonant effects. It requires 
an additional trajectory simulation (one without the resonance terms) 
than the indirect method previously used by the author [D 2 ] in analyzing 
2 k hour data. However, since it evaluates these effects directly, and 
without assuming a specific longitude gravity field, it should prove to 
be a valuable check on the previous method. 

For the 12 hour (one day) resonant orbit, the next relevant harmonic 
^nm next to least n satisfying (32) for r = ) is H44. 

The orbit averaged equation for the acceleration of the longitude of 
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the ascending equator crossing, due to H44 on a 12 hour orbit [(67), 


with yLi=.‘Xj ] gives: 


X 


r- 

at 


2 

rad,/sid.day • 


We consider a partical trajectory with only H44 acting as a perturbation. 
The initial conditions for the trajectory are otherwise the same as in 
the Hag studies. The H44 constants are: 10 /t^ ^ 2.2-4'^ 


The X" value 

H- 

is about 100 times a likely realistic value (see Figure 8) to provide 
an acceleration which will not be substantially "lost" in the roundoff 
error of the trajectory generator. 

From (DI2) the long term acceleration should be: 

X = -4.7S-3A»'f' 

Trt- ^ 2 

rad./sid.day , 


The following ascending equator crossing data is from a I6 day 
numerically integrated particle trajectory, in the presence of pertur- 
bations due to H44 only (whose constants are given above) . 


Time 

Longitude 


i 

e 

(Hrs.) 

(Begs.) 

(E.R.) 

(Degs. ) 


0. 

0. 

k . 1645 

30.0 

10"^ 

191 A82 

-0.08736 

4.1647 

30.0 

10"7 

382.975 

-0.3^907 

4.1648 

30.0 

10"'^ 

TABLE D5 

- Numerical Trajectory Data for 
Disturbed by H44 Only. 

a 12 Hour Earth 

Satellite 


(D12) 


(DI3) 
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In the above trajectory, the longitude is given approximately as: 


^ -Ti- ^-c-I71-7zI/s7.2S6 = 

f*.£Ay -»/ / 

rad./sid.day^ « (D 14 ) 

Comparison of (DI3) with (Dlh) again shows a very small discrepancy 
between the theory and the simulated data; of the order of 0.05^. 

16 HOUR ORBIT 

The 16 hour resonant orbit will undoubtedly be the most promising 
one for geodetic purposes, of the two synodic day orbits (see TABLE 5 )* 

For this orbit, the leading relevant harmonic is H43 (having the least 
n satisfying (32) for \~ ■= 3 ^ Ti s X. ). The orbit averaged equation for 
the ascending equator crossing due to H43 on a 16 ho\rr orbit [(67)7 
with ji-= 2 y p — I ] gives: 

radysid.day^ « (DI5). 


We consider a I6 hour orbit (a^ = 5 . 04^9 E.R.) with only H43 acting as 
a perturbation. The initial circular orbit has an inclination of 50 ^ 

A c 

— o . The H43 constants 

. o 

are: J43 = -2.2098x10“^, = <9 . Again^J43 is chosen about 100 

times realistic in order to overcome roundoff error in the comparison 



The following ascending equator crossing data is from a numerically 
integrated partical trajectory, in the presence of perturbations due to 
H43 only (with constants given above). 


TIME 

LONGITUDE 

a V 

i 

e 

(Hrs.) 

(Degs. ) 

■ (E.R.) 

(Degs.) 


0. 

0. 

5.0H9 

30.0 

10~® 

191.i^-7i^- 

0.02335 

5.0449 

30.0 

io“® 

382.9lj-5 

0.09107 

5.0448 

30.0 

CO 

1 

0 

1 — 1 


TABLE D6 - Numerical Trajectory Data for a I6 Hour Earth Satellite 
Disturbed by H43 Only. 


In the above trajectory the longitude acceleration (over the I6 
day period) is given approximately by (D 4 ), or: 




i^1L\00 

y 


rad./sid.day 


Again the discrepancy between measured (DI7) and theoretical (DI6) 
accelerations is about 0.15^. 

We note that for the weak resonances iX rad./sid.day ) 

the accuracy of the determination of the acceleration may be improved 
consideraily by utilizing all 2r equator crossings. In doing so^ of 
course^ care must be taken in choosing- the correct foirni of the condition 
equation (67); using^=r/ for descending and ^^ 2 . for ascending 
equator crossings. 

In summary^ the results of this section combined with the reported 
accuracies of the tracking of early bird^ LE>2J^ appear to make all the 
orbits in Table 3 amenable to simple geodetic analysis for resonant 


(D17) 
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longitude gravity harmonics in the earth's field. (See conclusions.) 
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APPENDIX E 


LIST OF SYMBOLS 


G 

X 


0 

i> 

s 


r 


i 


t 


& 


Location of the Greenwich Meridian 

Generally, geographic longitude of the ascending equator 
crossing of the satellite (also referred to as Xo , with 
reference to the beginning or near the beginning of the 
dynamics). Can also refer to the geographic longitude 
of the descending equator crossing, 
when subscripted, refers to orbit number 
Geographic latitude of the satellite. 

Position of the satellite ( % is its position at the 
start of the dynamics or at the reference ascending 
equator crossing) . 

Number of orbits in the synodic period of the resonant 
satellite (n* sidereal days); an integer. 

An integer (positive, negative or zero) the number of 
global circuits of the resonant orbit's ground track^ in 
its synodic period. When subscripted, alsOj integers of 
the rational fraction expressing the commensurability of 
the resonant orbit with the earth's rotation. 

Time from the beginning of the dynamics or some arbitrary 
zero. 

Argument of the ascending node, central angle in the orbit 
from the ascending node to the satellite (also subscripted 

^0 and ^ ) . 
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LIST OF SYMBOLS, COOT. 


Inertial longitude excursion of the satellite from its 
ascending node 

Inclination of the satellite’s orbit (also used as a diommy 
integer-subscript) . 

The earth’s inertial rotation rate. 

Geograpnic excursion of the satellite from an initial 
position 

Orbital revolution rate of the satellite. 

The inertial longitude turned by the earth since the 
start of the dynamics 

The distance from the center of the earth to the satellite. 
Perturbation forces in the radial, longitudinal and lati- 
tudinal directions, acting on the satellite 
The azimuth of the satellite’s trajectory. 

The synodic period (in sidereal days) of the resonant orbit 
satellite, an integer. The number of integral days for 
the stationary groimd track to begin to repeat itself. 

The semimajor axis of the satellite’s orbit^also sub- 
scripted is the semimajor axis in units of 

earth radii) . 

The tangential (along track) component of the perturbation 
force F. 

A perturbation force (when superscripted, a non-dimensional 
forc^ When subscripted nmp, an inclination force function . 
The total energy (kinetic and two-body gravitational potential) 


of a satellite. 


LIST OF SYMBOLS, COKT. 


n 

tr 

At 

c ) 

-r 

• « 

o 

% 


A small, but finite, change not sub or superscripted) 

The earth's Gavissian gravitational constant. 

The r orbit average of the quantity ( ), except in 

reference to -J” (see below). 

An increment of time^ in units of synodic days (not neces- 
sarily integral) 

An increment of time, in units of sidereal days. 

^ I } with the time increment in units of sidereal days. 

TT 

The period of the satellite's orbit ( ^To or is 

the period in sidereal days). 

vith the time increment in units of sidereal days, 
A nondimens ional semimajor axis change rate: — ^/o, ; 

m 

the units of arei ( sidereal days) 

The gravitational potential function of the earth. 

Signifying the gravitational harmonic teimi of order n and 


JL^c ) 
TP^ 




c 





power m. 

The mean equatorial radius of the earth. 

The associated Legendre polynomial of order n and degree m. 
The amplitude and phase of the non-normalized gravity 
harmonic H 

nm 

The normalized amplitude, cosine coefficient and sine 

coefficient of the gravity harmonic 

The ground track longitude span of half an orbit. 
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I 



I 



W «r 4 Vv/ 

P 

X' 

R-H 

« 

J 

F tn, F 

nmf> y Hinp 

u 

e-y^eC 3 
Y 

ic 

K»j 

c.c 

k' 

i 


When subscripted by integers and/or superscripted, 

a or b are constants 

Work done by perturbation forces. 

Case numbers 

A positive integer, or zero 
A dummy variable 

The mean anomaly of the satellite at the beginning of 
the dynamics , 

Keplerian elements; eccentricity (not used as a sub- 
script), argument of perigee, mean anomaly, true 
anomaly, and right ascension of the ascending node. 

The real part ^ ^ 

0 , except where used as a dummy integer. 

An inclination force function (when superscripted, the 
real value of the force function). 

A disturbing potential function 
s'- ^ , e = 2.718... 

The location of the Vernal Equinox. 

A positive integer or zero. 

Phase angles of the disturbing function. 

Phase and Frequency constants of the disturbing function. 
An integer or zero. 


[5] 

The number of orbits per day in Allan *s Resonance Theory 
As a subscript or in $ O j 

refers to a descending equator crossing^ fi= 2 or even 
refers to an ascending equator crossing. 
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, s 








‘nM 


Xcr Xe,5 ; X e^t 


' — /'y 
CA^ 
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Non-normalized cosine and sine gravity harmonic coefficients. 
Standard deviation (or estimate of it) 

Dominant , subdominant , 

An inclination dependent harmonic forcing function usually 
assumed constant over a libration period. 

An equilibrium longitude; stable^ and unstable. 

The minimum libration period of a resonant orbit. 

A frequency or angular rate. 

An inclination independent harmonic forcing function. 
Harmonic forcing functions. 
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